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SUMMARY 

A method is presented for measuring total electron content 
by observing the Faraday effect on satellite beacon transmitters 
at a station located on the magnetic equator. This method is ap- 
plied to data observed at Huancayo, Peru, between September 1961 
and February 1962. The total electron content at the magnetic 
equator was found to have nearly a 10 to 1 diurnal variation with 
a maximum of 4 x 10 17 electrons/m 2 near 1500 hours. Values of 
the thickness parameter (j Ndh)/ / N nia x for the observational period 
show an anomalous increase around layer sunrise which is at- 
tributed to a departure from thermal equilibrium. 
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SATELLITE OBSERVATIONS OF 
THE EQUATORIAL IONOSPHERE** 


by 

L. J. Blumle 

Goddard Space Flight Center 


INTRODUCTION 

Faraday rotation techniques have been used in conjunction with lunar radar investigations by 
other workers (References 1 and 2) to measure the total electron content of the ionosphere at mid- 
latitudes. The ambiguity in the number of polarization rotations was resolved by using two closely 
spaced frequencies. More recently beacon satellites have been used to determine total electron con- 
tent with the dispersive Doppler (also called differential Doppler) and Faraday rotation techniques 
(References 3 to 6). The total number of rotations has been estimated with reasonable accuracy by 
estimating the small number of rotations along a ray path that approaches transverse conditions 
(References 3 and 7). For a station located on the geomagnetic equator, there is a time during each 
satellite pass when the ray path and geomagnetic field are perpendicular and there are essentially no 
rotations. 

FARADAY ROTATION 

If it is assumed that both magneto-ionic modes follow a common ray path determined by the iso- 
tropic indices, the total rotation, in radians, is given by: 

fi = -X J K ->-) ds - (1) 

Assume an index of refraction of the form (Reference 8): 

jj? ~ \ - i + y l * (2) 


♦ This worlT was don e while the author was at the Ionosphere Research Laboratory, Pennsylvania State University, and was supported by the 
National Aeronautics and Space Administration under Grant NSG 114-61. 

+ A similar paper was published in J. Geopbys. Res . 67(12):4601-4605, November 1962. 
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where 



in which 

f = wave frequency, 
f N = plasma frequency, and 

f L = the component of the gyro frequency along the ray path. 

A binomial expansion of Equation 2 neglecting all but first order terms yield the result 


n 



RAY 

PATH 


XY, 


V 


ds . 


(3) 


For a spherically stratified ionosphere, Equation 3 becomes, upon substitution for x and Y L : 


n 


2.36 X 10 4 


B l secy 



f 2 


(4) 


where 


x 

N 

h 


B l sec y 


component of the magnetic field along the ray path in webers/m 2 , 
zenith angle of the ray path, 
electron density in electrons/m 3 , 
height, and 

a weighted mean for the range of integration. 


When the second-order terms of the binomial expansion of Equation 2 and refraction (Reference 
9) are considered, the correct second-order expression is: 


n = 


Xo xy l 1 + T 


RAY 

PATH 


x x tan 2 y 0 


dh , 


( 5 ) 


where x 0 is the zenith angle of the ray at the ground (Reference 5). 


2 



At 54 Me the refraction correction may be neglected for small zenith angles (x 0 < 50 degrees). 
Furthermore, the integration may be carried out along a straight-line ray path without appreciable 
loss of accuracy. If a reasonable distribution of ionization is assumed and refraction is neglected, 
Equation 5 becomes 


n 


2.36 x 10 4 


B l sec N dh 



- 

(f 0 F 2 \ *1 

f2 

1 - 0.35 1 

If) J 


( 6 ) 


where f 0 F2 is the critical frequency of the F region of the ionosphere. 

It is important to bear in mind the assumptions made in the derivation of Equations 3 through 6: 

1. The propagation path satisfies the quasi- longitudinal (QL) approximation of the Appleton- 
Hartree formula (Reference 8). 

2. In the region where the index of refraction differs from unity, the wave frequency is much 
greater than the electron gyro-frequency, the plasma frequency, and the electron collision 
frequency. 

The height at which B L sec~7 is the weighted mean for the range of integration can be determined 
if the N-h profile is known. A study of published N-h profiles indicates that the mean height is 400 ± 25 km 
for integration to 1000 km (Reference 10). For the results reported here, B L sec x has been evaluated 
by using a spherical harmonic analysis of the main field. The coefficients used in this analysis were 
obtained from Reference 11. 


TRANSVERSE PROPAGATION REGION 

It is important to consider the behavior of the polarization through the quasi-transverse (QT) re- 
gion. Figure 1 illustrates the mode behavior in the QT region and is calculated for a wave frequency 
of 54 Me with a typical daytime ionosphere model. Under QL conditions, the two magneto-ionic modes 
are circularly polarized and the phase path difference between modes, expressed in wavelengths, is 
the total number of half- rotations of the plane of polarization. 

As the transverse region is approached, the mode polarizations become elliptical and their re- 
spective indices depart from the QL approximation. The QT region may then be identified from the 
resulting elliptical wave polarization. Because of this resulting elliptical wave polarization, the 
phase path difference cannot be simply interpreted as rotation of the plane of polarization (Reference 
12). The typical time scale for the 54-Mc ray path from an earth satellite in a 1000-km circular orbit 
to sweep through this transverse region indicates that the QL approximation (Equation 6) may be ap 
plied accurately at 54 Me to deduce total rotation unambiguously throughout a satellite pass except 
for about 6 seconds on each side of the transverse point. 
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OBSERVATIONS 

The satellite. Transit 4A (1961 Ol), was 
launched into an approximately circular orbit 
inclined 67 degrees to the equator. The satellite 
is oriented along the local geomagnetic field. 
The 54- Me radiation pattern is such that linear 
polarization with stable orientation is radiated 
near the magnetic equator. 

Since September 1960 this satellite has been 
observed by a station located on the magnetic 
equator near Huancayo, Peru (Reference 13). 
The signals are received on fixed folded-dipole 
antennas by fixed- gain crystal controlled re- 
ceivers, the outputs of which are recorded on 
magnetic tape along with timing marks from a 
clock. 

Figure 2 is a pen record produced from the 
magnetic tape of a typical midday pass at 
Huancayo during local summer. The top chan- 
nel displays the timing marks, and the other two 
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Figure 1 —Parameters for 54-Mc propagation in a 
typical daytime ionosphere. 
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channels are linear displays of the amplitude received on two oriented dipole antennas. The spacing 
between nulls on each channel corresponds to a rotation of the plane of polarization of v radians; for 
the northbound pass shown here the rotation is from north into west. 

At 1854 UT, the satellite is in the southern magnetic hemisphere and the propagation is QL. Near 
1854:45 UT, the wave polarization is elliptical, and this is identified as the transverse region QT. This 
elliptical polarization is obseiwed for several seconds during each daytime pass, as Figure 1 indicates. 


TOTAL ELECTRON CONTENT 

Total electron content is measured by simply counting the number of rotations observed between 
the transverse point QT and some later time when the satellite is in the QL region. The satellite lo- 
cation at this time is obtained from the orbital ephemeris and used to compute the term B L sec * enter- 
ing in Equation 6. From Equation 6 it is possible to obtain values of total electron content with better 
than 10 percent accuracy. 

The values obtained from several applications of this technique to each satellite pass are then 
averaged (standard deviation less than 5 percent) and plotted as a function of local mean time (Figure 3). 
The nodal precession of this satellite results in passes over the station occurring, on the average, 
earlier each day so the diurnal variation in total electron content could be measured over a six-week 
interval. The data plotted in Figure 3, observed from September 1960through January 1962, correspond 
to several cycles of quiet- day observations. 


Figure 3 — Diurnal variation of total elec- 
tron content at the magnetic equator, Sep- 
tember 1961 -January 1962. 



LOCAL MEAN TIME 
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No seasonal variation was evident on the data and this is not unusual when we consider that the 
solar noon zenith angle was less than 12 degrees over this observational period. 

The nearly 10 to 1 diurnal variation is the outstanding feature of the values of total electron con- 
tent. A minimum of about 0.5 x 10 17 electrons/m 2 at 0300 hours and a maximum greater than 4 x 10 17 
electrons/m 2 near 1500 hours may be contrasted with mid-latitude values for the same period which 
show a broad maximum of about 2 x 10 17 electrons/m 2 near 1600 hours and nearly a 2 to 1 diurnal 
variation (W. J. Ross, private communication). 


THICKNESS PARAMETER 

Since the major contribution to the total electron content is in the F region, the thickness param- 
eter £fNHh )/N max is closely associated with the scale height of the F region. If the N-h profile is a 
Chapman function with constant scale height, the total content is given by 


where 


Ndh 


4.13N H 

max 


( T . +T i) 

mg 


( 7 ) 


( 8 ) 


in which 


K = Boltzmann's constant, 1.38 x 10" 23 joule/°K, 

T e = electron temperature 
T. = ion temperature 
m = molecular mass, and 
g = acceleration of gravity. 

For N-h profiles other than that of the familiar Chapman expression, Equation 7 will be modified 
by a different constant. 

The thickness parameter (j Ndh) / /N max is plotted as a function of local mean time in Figure 4 along 
with theoretical values of neutral gas temperature based on the value of the 10.7-cm solar flux appro- 
priate to the observational period (Reference 14). Values of N majt were obtained from extrapolation 
in solar time of vertical incidence values observed at Huancayo (Instituto Geofisico del Peru, private 
correspondence). The hourly mean curve of the thickness parameter shows a distinct anomaly near 
sunrise. Measurements at Lima, Peru, of the effective scattering cross-section of free electrons by 
the radar backscatter technique (Reference 15) indicate that only during this sunrise period are de- 
partures from thermal equilibrium significant (T e /T. = 2) . It is suggested that the anomaly in the 
thickness parameter near sunrise is due to this departure from thermal equilibrium. Bauer 
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Figure 4— Diurnal variation of the equatorial thickness parameter [JN dh)/N 
September 1961 -January 1962. ' ' 


(Reference 16) and Evans and Taylor (Reference 2) have observed similar departures in the effective 
scale height during the sunrise period. 
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